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Abstract 



We consider the two-loop four-point amplitude in J\f = 2 super QED, and show that there exists 
an approximate recursive structure similar to that captured by the ABDK/BDS ansatz for MHV 
amplitudes in TV = 4 super Yang-Mills. Furthermore, we present a simple relation between the box 
coefficients of one-loop photon MHV amplitudes in (super) QED, and sums of box coefficients of 
one-loop MHV amplitudes in (super) Yang-Mills. 
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1 Introduction 



One of the important realisations of the past decades is that physical observables in quantum gauge 
theories are far simpler than one would expect from Feynman diagrams. For instance, the Parke- 
Taylor formula [1] for the maximally helicity violating (MHV) scattering amplitudes in colour-ordered 
Yang-Mills theory at tree level resums large numbers of Feynman diagrams into a stunning one-line 
expression. Such intriguing simplicity persists at the quantum level, culminating perhaps in the higher- 
loop iterative structures discovered in the loop expansion of maximally supersymmetric Yang-Mills 
(SYM) in [2,3]. 

The perturbative expansion of supergravity theories is also full of surprises. At tree level, there 
are interesting relations between amplitudes in Yang-Mills and in gravity, starting with the KLT 
relations [4] and continuing with the recent solution of the BCF recursion relations [5, 6] for general 
relativity [7,8] found in [9], which expresses amplitudes in maximal supergravity in terms of sums 
of squares of AA = 4 SYM amplitudes. Both KLT formulae and the relations of [9] have echoes in 
the expressions for the one-loop box coefficients [10-13]. Most importantly, there is now mounting 
evidence of the remarkable similarities between M = 4 SYM and N = 8 supergravity, leading to 
the conjecture that the J\f = 8 theory could be ultraviolet finite, which is supported by multi-loop 
perturbative calculations [14-17]. 

In describing the remarkable web of regularities and similarities between the perturbative expan- 
sions of gauge theory and gravity, Quantum Electrodynamics (QED) has its own place in the story. For 
example, multi-photon amplitudes in QED (with at least eight photons) have in common with max- 
imally supersymmetric Yang-Mills and supergravity the no-triangle (and no-bubble) property. This 
is the statement that all one-loop amplitudes can be written as sums of box functions times rational 
coefficient^- This property was proven for Af = 4 SYM in [18], conjectured for Af = 8 supergravity 
in [10, 11, 19,20] and subsequently proved in [21,22]. Recently it was found in [23] that a similar 
statement holds for photon amplitudes in QED. We also mention the interesting connections found 
in [16] and [21,23] between the unexpected cancellations in one-loop scattering amplitudes, and the 
large- z behaviour of tree amplitudes observed in [7,8,22,24,25]. In unordered theories such as gravity 
and QED these cancellations are amplified by the summation over different orderings of the external 
particles. 

Two more interesting facts are worth mentioning. Firstly, the one-loop MHV and four-point two- 
loop photon amplitudes in M = 2 SQED have a uniform degree of transcendentality, i.e. at one and 
two loops, only terms with total polylogarithmic weight equal to 2 and 4 appear, respectively [26]. A 
similar fact has been recently observed in [27, 28] in the one- and two-loop graviton MHV amplitudes 
in maximal supergravity. Furthermore, the M = 2 SQED result for these amplitudes can be obtained 
from the corresponding N = 1 SQED result by keeping only terms with maximal transcendentality 
(and no ratio of kinematical scales), leading to the speculation that maximal transcendentality [29] 
could be a feature of all maximally supersymmetric theories. Moreover, slightly departing from the 
realm of scattering amplitudes, we would also like to recall the somewhat puzzling "simplicity" of the 
three-loop electron anomalous dimension [30]. Here, numerically large cancellations occur between 

1 One-loop photon amplitudes in (S)QED are somewhat special since they are both infrared and ultraviolet finite. 
This implies particular relations between the box coefficients, since the infrared divergences must cancel. 
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different diagrams, a fact which is due to the breaking of gauge symmetry at the diagrammatic level, 
see [31] for a prescient and enjoyable discussion of this point. 



It is therefore natural to ask to which extent the simplicity found in the perturbative expansion 
of amplitudes in supersymmetric Yang-Mills and supergravity persists in (S)QED. We are fortunate 
to have a large number of analytic amplitudes at our disposal to test this. The one-loop four-photon 
amplitudes for massless and massive fermions were first computed in [32,33]. Corrections to light-by- 
light scattering at two loops were determined about fifty years later in [34] using the modern unitarity 
method [18,35]. The four-point results of [34] were confirmed in [26] and extended to N = 1 and 
N = 2 SQED by analysing the tensorial structure of the amplitudes found in [32,33]. In [36], analytic 
expressions for one-loop MHV photon amplitudes for an arbitrary number of photons were calculated 
with the help of the off-shell currents found in [37] . In [38] , analytical results for all six-photon QED 
amplitudes were given whilst in [39], formulae for n-point MHV amplitudes in QED, scalar QED and 
M = 1 SQED were obtained, together with the analytical results for the six-point NMHV QED and 
M = 1 SQED amplitudes, which confirmed earlier work of [36,38]. 



At tree level, the simplest nonvanishing scattering amplitude one encounters in massless QED is 
the MHV amplitude with n photons and two fermions H 

Amv(^,i + ,2+,...,r,...,n+) = i^nT||- (1.1) 

{qq) 2 fj[{ql){iq) 

. (qi) 3 (qi) sr^ (qq) 

where the fermion q and the i th photon have negative helicity, and all the other particles have positive 
helicity. Equation (jl.l|) shows two important features. Firstly, the MHV amplitude in QED is given 
by a compact, one-line expression, see the first line of (11.11) . Furthermore, this amplitude can be 
derived by summing over permutations of colour-ordered amplitudes in Yang- Mills where the photons 
are replaced by gluons with the same helicitiesH This is explicitly shown in the second line of (II. ip . 
where each term in the sum over permutations V{1, 2, . . . , n} is equal to a colour-ordered Yang-Mills 
MHV amplitude with n gluons and two fundamental fermions q and q. 



This observation leads directly to the first result we present in this paper. We will discuss how 
the one-loop MHV amplitude of photons in supersymmetric and in pure QED can be derived directly 
by summing over appropriate permutations of the corresponding result for gluon MHV amplitudes in 
supersymmetric or pure Yang-Mills theory. As we mentioned before, one-loop photon amplitudes in 
(S)QED can be written in terms of (the finite parts of) box functions for n > 8. We will therefore 
show that the box coefficients of the Yang-Mills amplitudes, summed over appropriate permutations 
of the external gluons, directly give the box coefficients of the QED amplitudes. We will also outline 
the proof of this interesting fact, based on a MHV diagram calculation [40-44]. 



The second observation we make in this paper is aimed at uncovering possible cross-order relations 
in the perturbative expansion of N = 2 SQED. The first example of iterative structures was found in 
planar N = 4 SYM for the four-point MHV amplitudes in [2]. In a subsequent paper [3], Bern, Dixon 

2 There is no tree-level photon amplitude corresponding to the gluon MHV amplitude in Yang-Mills. 
3 See [51] for a discussion of this important feature. 
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and Smirnov (BDS) put forward a conjecture for the all-loop resummation of the planar n-point MHV 
amplitudes in M = 4 SYM, which has been tested up to three loops in the four-point case [3] and 
up to two loops in the five-point case [45,46]. However, in [47] it was realised that the BDS ansatz 
is incomplete at least for a large number of external gluons. Specifically, a direct calculation of the 
two-loop six-gluon MHV amplitude in M = 4 SYM performed in [48] showed that the BDS ansatz 
breaks down, and has to be amended by adding a dual conformal invariant remainder function [49,50]. 

Motivated by this, we will consider the four-photon MHV amplitude at one and two loops in 
maximally supersymmetric N = 2 SQED, and test the possibility that the two-loop amplitude could be 
written as a polynomial in the one-loop amplitude. One important difference compared to Yang-Mills 
is that in QED the one- and two-loop four-photon amplitudes are finite. Thus, one lacks the guiding 
principle of the exponentiation of infrared divergences, which is central to the all- loop ABDK/BDS 
ansatz. Despite this, we find that, quite surprisingly, the real part of the two-loop four-photon MHV 
amplitude in maximally supersymmetric QED is not exactly given but well approximated (in a wide 
kinematic region) by a polynomial in the one-loop MHV four-photon amplitude. We also discuss the 
limitations of such an approximate formula. 

The rest of the paper is organised as follows. In Section 2 we present the relationship mentioned 
earlier between the box coefficients of one-loop MHV (S)QED amplitudes and sums of permutations 
of box coefficients of the same amplitudes in (S)YM, and prove it using one-loop MHV diagrams. In 
Section 3, after reviewing salient features of the BDS ansatz, we investigate approximate recursive 
structures for MHV four-photon amplitudes in J\f = 2 SQED. 



2 One-loop photon amplitudes in massless (S)QED 

In this section we wish to comment on a simple relation between massless scalar QED and pure Yang- 
Mills amplitudes, as well as a similar one between M = 1 SQED and M = 1 SYM amplitudes. Similar 
relations, based on certain permutation sums of gluon amplitudes, are known for tree amplitudes [51] 
and one-loop amplitudes (see e.g. [52,53]). 

We start by considering the expressions for the scalar QED and M = 1 SQED photon MHV 
amplitudes at one loop. These amplitudes were computed in [39] , and are given by 





(2.1) 



+ i 
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where 



lar o (13) (14) (2 3) (2 4) " (34)"- 4 

(3 4)2 11 {3im) • 



Let us explain the notation employed in (|2.ip . Firstly, the sums in (|2.ip are over permutations V of 
the massless states inside the curly brackets. Secondly, the function B 2me appearing in (|2.ip is the 
finite part of the two- mass easy scalar box functions F [42,54], 

F 2me (s,t,P,Q) = -±[(-s)-* + (-t)-*-(-P 2 )-*-(-Q 2 )- £ ]+B 2me {s,t,P,Q), 

where 

B 2me (s,t,P 2 ,Q 2 ) = Li 2 (l -aP 2 ) +Li 2 (l- aQ 2 ) -Li 2 (l- as) -Li 2 (l -at) , 

with 

P 2 + Q 2 -s-t 
a: = P 2 Q 2 -st ■ (2 - 4) 

As usual, s := (P + p) 2 , t := (P + q) 2 , with p + g + P + Q = 0, where p and (/ are the massless 
legs (sitting at opposite corners, in the two-mass easy boxes), and P and Q the massive legs. The 
arguments of the box functions appearing in (|2.ip are the kinematical invariants Sj...,- := (&« + • • - + kj) 2 . 



In Figure 1 we provide a representation of the box function appearing in the second line of (|2.ip . 
The massless legs correspond to the positive helicity photons 3 + and 4 + . The negative helicity photons 
1~ and 2~ are always part of (different) massive corners P and Q, which contain m — 3 and n — m + 
1 legs respectively. The combinatorial coefficients appearing in (12. ip correspond to the number of 
permutations of the positive helicity photons inside P and Q (which obviously leave the box function 
invariant). 

In (|2.2|) and (|2.3|) , we have multiplied the result of [39] for rf scalar by a factor of 2 to account for the 
fact that we are working with complex scalar fields. Finally, let us stress that the amplitudes given in 
(|2.ip are infrared and ultraviolet finite. Because of Furry's theorem, they are nonvanishing only for n 
even. 



We now turn to the corresponding planar MHV amplitudes in Yang-Mills. In J\f = 4 SYM they 
were first derived by Bern, Dixon, Dunbar and Kosower in [18] using unitarity [35] and collinear limits, 
and later confirmed in [42] using one-loop MHV diagrams. They have the following form 

i-l j-i 

A*= A = AT £ £ F 2me (p,q,P,Q), (2.5) 

p=j+l i?=i+l 



and 



^■(i + .....r,-,r...... + ) = -« (ia>(2 ^.. (Bl) , (2.6) 
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Figure 1: The two-mass easy box function appearing in (|2.1f) . The momenta p = k% and q = are 
null, whereas P := k\ + ks + • ■ ■ + k m and Q := hi + k m+ i + • • • + k n are massive. The one-mass 
box function in (12. ip is obtained by setting m = n, so that the top right corner becomes massless (and 
contains only the momentum foj. 



is the tree- level amplitude, given by the Parke- Taylor formula [1]. 

The one-loop MHV amplitude in J\f = 1 SYM was presented in [35] and rederived in [43] using 
MHV diagrams. The contribution to the amplitude of an N = 1 chiral multiplet running in the loop 
is given by the following compact formula, 

<= 1 (l+,...,r,... ) r ! ---,n+) = ^ l< X=l } l P J q B 2 ™{p, q ,P,Q) + ... , (2.7) 

p=j+l q=i+l 



where 



and 



i^ =1 ]% = \ATb%, (2i 



,{i p){j q){iq) {jp) 
(ij) 2 (pq) 

The dots in (|2.7p stand for triangle and bubble functions, which do not enter our discussion@ 



Lastly, the one-loop n-point non-super symmetric Yang-Mills MHV amplitudes were computed 
in [35,44,55] and confirmed in [56] using generalised unitarity [57,58], with the result 

-i j'-i 

AT l ™(l + ,---,r,...,r,...,n + ) = E [c SCalar ]%B 2me (p,q,P,Q)+--- , (2.10) 

p=j+l q=i+l 



4 This is because of the no-triangle property of QED amplitudes [23], which ensures that the bubble and triangle 
coefficients vanish. 
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where 

[c scala %{ = \ATm 2 - (2.11) 

The dots in (|2.1Up stand for triangle and bubble functions, as well as for the rational terms of the 
amplitude, which are not relevant in the following. 

Next we wish to expose a simple relation between the coefficients cf caiar/.A/=i Q f the two-mass easy 
box functions in the expression for the MHV photon scattering amplitudes in (S)QED, given in (|2.ip . 
and the corresponding coefficients of the same box function, c scalar Mf=i^ Q f the MHV gluon amplitudes 
in N = 1 and pure Yang-Mills in (|2.8I) and (|2.1ip . In order to match a gluon amplitude to a target 
(S)QED amplitude, we introduce appropriate sums over permutations as follows. With reference to 
the box function in Figure 1, we hold the two massless legs fixed, and sum over permutations of the 
gluons appearing at the massive corners. The sum can be performed using the eikonal identity [51], 
see also (11.11) . to get 



EE 



<12) 4 (12) 4 (34)- 4 



j (3 2) (2 (m + 1)> ■ ■ ■ ((n - 1) n) (n 4) (41) (15) • • • (m3) (3 1) (1 4) • • • (3 n) (n 4) ' 

where V\ := V{2, m + 1, . . . , n} and V2 '■= V{1, 5, . . . , m} are permutations of the massless legs in the 
massive corners of the box function. Multiplying (|2.12|) by 634 we recover the expression for [c^ =1 ]34 
given in A2.2H . A similar argument runs for the one- mass box coefficients [d scalar ]p g , with the only 
difference that the sum in (|2.12p is over one set of permutations rather than two. 



One can arrive at the same conclusion by performing a one-loop MHV diagram calculation akin 
to [42-44] with MHV rules adapted to QED as done at tree level in [41]. For a one-loop MHV photon 
amplitude we have to glue two tree-level MHV vertices with two internal scalar propagators, and 
perform an appropriate loop integration [42]. We will not give details of the calculation because we 
can recycle results from [42-44]. The crucial observation is that the only diagrams contributing are 
those where all gluons are external and the two internal legs of each MHV vertex are either scalars 
or fermions with opposite helicity. This also implies that the two external negative helicity gluons 
must belong to different MHV vertices. The relevant MHV QED tree amplitudes can be obtained 
from the corresponding MHV tree amplitude in QCD with n — 1 positive helicity gluons, one negative 
helicity gluon and two fermions (scalars) of opposite helicity and summing over the n! permutations 
of the n gluons. Writing the QED MHV vertices with two fermions and n gluons in terms of QCD 
tree MHV vertices, see the second line of (jl.ip . reduces the calculation to sums of permutations of 
MHV one-loop diagrams for MHV amplitudes in pure Yang-Mills and N = 1 SYM [43,44]. It can 
be easily seen that this reproduces exactly the observations made earlier in this section on the box 
coefficients of the one-loop (S)QED MHV amplitudes. Triangle and box coefficients are guaranteed 
to vanish because of the no-triangle property [23]. 



Finally, we observe that this has implications for the twistor-space localisation properties of the 
coefficients, which are inherited from those of the (S)YM amplitudes, i.e. the coefficients localise on 
sets of two, possibly intersecting lines in twistor space. It would be interesting to see whether similar 
structures appear in non-MHV amplitudes and at higher loops. 
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3 Approximate iterative structures in Af = 2 SQED 



Motivated by the existence of iterative structures for amplitudes in M = 4 SYM, we have investigated 
the possible existence of recursive-like structures for MHV amplitudes in the maximally supersymmet- 
ric M = 2 SQED theory! Before discussing our results, let us briefly review the iterative relations in 
M = 4 SYM [2,3]. It was shown in [2] that the two-loop four-point MHV amplitude in M = 4 SYM 
satisfies an intriguing cross-order relation, 

M?\e) - \{M { l\e)f = fWie)MP(2e) + C<a+0(e). (3.1) 

Here is the helicity-blind function obtained by taking the ratio between the L-loop MHV 

amplitude and the corresponding tree amplitude. Furthermore f^ 2 \e) = — (C2 + C3 e + C4 e2 )? an d 
C( 2 ) = -(5/4)C 4 . 



In [3], a resummed, exponentiated expression for the scalar function M. n was proposed, and checked 
explicitly in a three-loop calculation in the four-point case. The BDS conjecture is expressed as [3] 

00 00 
M n := 1 + J> L M L) (e) = ^ V [Y j a L {f^{e)M^{Le) + C^+E^{e))] , (3.2) 

L=l L=l 

where a = [g 2 iV/(87r 2 )](47re~ 7 ) € . Here f^ L \e) is a set of functions, 

f (L) (e) := ti L) + f[ L) e + A L) e\ (3.3) 

one at each loop order, which appear in the exponentiated all-loop expression for the infrared diver- 
gences in generic amplitudes in dimensional regularisation [59] (and generalise the function in 
(13. ip ). In particular, fff^ = 7^ /4, where is the cusp anomalous dimension, related to the anoma- 



lous dimension of twist-two operators at large spin. Importantly, the constants C^ L \ Jq L \ f[ and 
on the right hand side of (|3.2p do not depend either on kinematics or on the number of particles 
n. On the other hand, the non-iterating contributions depend explicitly on n, but vanish as 
e 0. 



BDS also suggested a resummed expression for the appropriately defined finite part of the n-point 
MHV amplitude, 

Tn = e F *° S , (3.4) 

where 

^n BDS (a)= \ lK {a)F^{U) + C{a) . (3.5) 

Notice that the entire dependence on kinematics of the BDS ansatz enters through the finite part of 
the one-loop box function, Fn (0). 

A similar analysis has been performed in [27, 28] for the four-point MHV amplitude in N = 8 supergravity, and 
highlighted a remarkably simple structure for the two-loop term in the expansion of logarithm of the helicity-blind ratio 
M M = 8 /Mtree- The functions appearing in the ratio were also found to have uniform transcendentality. 
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In analogy with the BDS ansatz, we would like to investigate the existence of cross-order relations in 
the four-point amplitude A^l - , 2~, 3 + , 4 + ) in J\f = 2 SQED. To this end, we consider a decomposition 
of the two-loop term in the expansion of this amplitude as 

[Mf ]ansatz = & [M^] ' + ^i'' + * , (3.6) 

where is the four-point one-loop MHV amplitude, and b, c and d have to be determined. 



The expressions for photon-photon scattering amplitudes at one and two loops entering (I3.6P are 
taken from [26]. The one-loop four-photon MHV amplitude in N = 1 SQED is given by 



M 



(i) 



-A[(X - Yf + 7T 2 



where 



X 



log ( — 

s 



Y 



log ( — 

s 



(3.7) 



(3.8) 



A few comments are in order here. Firstly, we notice that in the physical region s > and t, u < the 
expression (13. 7h is real. Outside this region, an analytic continuation is needed as the u- and t-channels 
develop a discontinuity. Secondly, we observe in (|3.7|) that ratios of kinematic scales such as t/s only 
appear as arguments of logarithms. Furthermore, all the functions appearing in the expression for 
have uniform degree of transcendentality equal to 2. 



The two-loop expression for the four-point MHV N = 2 SQED amplitude is still rather compact 
and simple. It is given by [26] 



M 



(2) 



-16Li 4 (y) + 8 YLi 3 (x) + 8V Li 3 (y) + 

45 



2 , 
-XYtt 2 
3 



+ m 



16Li3(x 



-Y 3 (Y- 
3 V 

' 3 



AX) 



-Y 2 (Y -3X) 



| it i| , 



where X and Y are defined in (j3.8l) and 

x := —t/s 



y ■-- 



-u/s 



(3.9) 



(3.10) 



As in (|3.7p . also in (|3.9|) there are no terms proportional to ratios of Mandelstam variables, and in 
(13. 9p we only have functions with transcendentality equal to 4. Therefore, we expect the coefficients 
b, c and d to have transcendentality and 2 and 4, respectively. In contrast to the one-loop amplitude 
the two- loop amplitude develops an imaginary part. 



The real part of (|3.9p can be recast in the following suggestive form, 



Re 



-16Li 4 (y) - 16Li 4 (x) + 8{X + Y) (hi 3 (x) + Li 3 (y)) + 4X 2 Y 2 - -XYir 2 



1 

24 L 



M 



(i) 



7T 



45 



-7T 



(3.11) 



(i) 



In (13. lip , we have re- written the result of [26] in a way that already incorporates the functions M 
and [.M^] 2 appearing in (|3.6|) . see the last line of (I3.11|) . In the following we will focus only on the 

(2) 

real part of M\ . 
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In order to find a set of coefficients b, c and d which best fit our proposed ansatz fl3.6f) . we build a 
system of three equations for some three random values of y and solve for b, c and d. We then plug 
the values of these coefficients in (I3.6P and match the ansatz against the real part of (13, 9p . In Figure 
2 the real part of f|3.9|) and our ansatz are plotted. The values of b, c and d used in the plot are given 
by the set of values / for y shown in Table 1 below, 



Coefficients 


I 


II 


III 


IV 


V 


VI 


VII 


b 


-0.0386 


-0.038 


-0.0387 


-0.0417 


-0.0412 


-0.0415 


-0.0412 


c 


-2.567 


-2.571 


-2.574 


-2.812 


-2.894 


-3.009 


-2.877 


d 


-5.784 


-5.894 


-5.938 


-11.46 


-14.689 


-25.098 


-14.489 


F{b,c, d)*= 2 


97.179 


187.129 


88.452 


49.555 


1.065 


25.554 


0.500 



Table 1: Values of b, c and d for different sets of values for y. The sets I -1 1 1 include points away 
from the boundary y = and y = 1, and the resulting coefficients b, c, d vary slowly from one set to 
another, unlike the case of the sets IV -VI, which include points near the boundary in y space. The 
set VII is obtained using the least square method. 

where the sets I-VI are| I = {0.3,0.4,0.5}, II = {0.272,0.342,0.482}, III = {0.25,0.35,0.45}, 
IV = {0.00001,0.10001,0.482}, V = {0.0006,0.0023,0.006}, VI = {0.0001,0.0003,0.0235}, while the 
values of (6, c, d) in column VII are obtained applying the least square method. We have carried out 
a similar analysis for the four-point MHV amplitude in N = 1 SQED, the expression of which can 
also be found in [26], and we report the results in Table 2 below, 



Coefficients 


I 


II 


III 


IV 


V 


VI 


VII 


b 


-0.0356 


-0.0356 


-0.0356 


-0.0277 


-0.027 


-0.026 


-0.0278 


c 


-1.829 


-1.828 


-1.827 


-1.062 


-0.113 


0.211 


-0.719 


d 


-79.04 


-78.99 


-78.98 


-60.99 


18.39 


32.54 


-44.57 


F(b, c, d) N=1 


756.0 


749.7 


746.5 


99.4 


1265.9 


1365.6 


24.1 



Table 2: Values ofb, c and d for different sets of values for y. The sets of y used are the same as in 
the TV = 2 theory shown in Table 1. 



Introducing 

F(b,c,d) := J dy(M?\y) - [A^f ]ansatz(y)) 2 , (3.12) 
and minimising F(b, c, d) gives the "best" values for b, c and d over all of the phase space. 

Let us summarise the outcome of this analysis. 

1. Unlike the case of the BDS ansatz for M = 4 SYM, we find that the coefficients b, c and d are 
not (kinematic-independent) constants, but have different values for different kinematic points. 

2. Having derived values of these three coefficients, we plot the M = 2 two- loop amplitude as well 
G Sincc (|3.9p is symmetric under x — » 1 — x, we only choose values for y from half of the phase space. 
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Figure 2: In this Figure we plot (the real part of) the right-hand side of (13. 9p . representing the four- 
point two-loop MHV amplitude in M = 2 SQED, together with our ansatz (|3.6|) . with y = —u/s given 
by set I. The two overlapping functions are indistinguishable in this plot. 

as our ansatz as functions of the ratio y. A particular example is presented in Figure 2, which shows 
a very surprising overlap of the two functions. 

3. In order to study more closely the functional forms of the two-loop amplitude and of our ansatz, 
in Figure 3 we present a plot of the difference between (|3.6p and (|3.1ip . which we could consider as 
the "remainder function" for this amplitudelll 



Tl 4 (y) := K<b\m¥ } J - lb [M ( 4>] z + cM\ L} + dj . (3.13) 

From Figure 3, we see that the two functions agree for a wide range of values of y, however the 
difference function has spikes as y — ► or y — ► 1. In these limits the best fit would be given by the 
second line of (|3.1ip in which case the disagreement would be proportional to a single power of a 
logarithm in x or y. The presence of this divergent behaviour at the extrema of the phase space shows 
that our ansatz is certainly incomplete. However, we find it quite remarkable that large deviations 
only appear at the boundary of the phase space. 

4. We have computed (|3.12p for both AT = 1 and M = 2 SQED and found the values 

F(b, c, d)^= 2 = 0.5 , F(b,c,d) u=1 =24.1. (3.14) 



This shows that (|3.6p gives the most accurate approximation of the real part of the two-loop amplitude 
in the case of M = 2 SQED. 

5. We also observe that the four-point MHV amplitude in M = 2 SQED can be derived from the 
M = 1 SQED amplitude by keeping maximally trascendental terms and deleting contributions which 
multiply ratios of the kinematics invariants. 



7 Notice that this function depends however on the choice of b, c and d we make in the ansatz (|3.6 
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Figure 3: This figure shows a plot of the M = 2 SQED remainder function, constructed as the difference 
between the two-loop MHV amplitude in N = 2 SQED and (|3.6|) using set I. 



In summary, the approximate iterative structures we have explored in this section are certainly 
not on the same footing as those in M = 4 SYM found in [2,3]. Nevertheless, we find it intriguing that 
part of the two-loop four-photon MHV amplitude (|3.9|) is captured by a polynomial in the one-loop 
amplitude (|3,7|) . It would be interesting if one could find exact iterative structures, written in terms 
of an appropriate QED remainder function, and explain them in terms of some underlying symmetry 
of the theory. It would also be interesting to find a Wilson loop interpretation of MHV scattering 
amplitudes in SQED, similarly to that found in J\f = 4 SYM [60-62]. 
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